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ABSTRACT 

We investigate the relationship between accretion rates and the spectral energy distributions (SEDs) of BL Lac 
objects, using a sample of objects for which published information on the host galaxies, emission-line luminosities, 
and peak frequencies and luminosities of their SEDs are available. The sample is composed of 43 BL Lac objects 
which have a relatively continuous distribution of peak frequencies. Under the assumption that the observed 
emission lines are photoionized by the central accretion disk, we use the line luminosities to estimate the accretion 
luminosities and hence accretion rates. We find that low frequency-peaked BL Lac objects (LBLs) span a wide 
range of accretion rates, whereas high frequency-peaked BL Lac objects (HBLs) cover a more restricted range of 
lower values. There appears to be a continuous distribution of accretion rates between the two subclasses of BL 
Lac objects. We find that the peak frequency of the SED, u p k, correlates with the accretion rate, approximately 
with the form oc A" 3 in HBLs and j/ p k oc A" 25 in LBLs, where A = Ly mes /c 2 . The peak luminosity of the SED 
is also correlated with A. These results suggest that the accretion rate influences the shape of the SED in BL Lac 
objects. They also support models which couple the jet and the accretion disk. We present a physical scenario to 
account for the empirical trends. 

Subject headings: accretion, accretion disks — black hole physics — BL Lacertae objects: general — galaxies: 
active — galaxies: jets — galaxies: nuclei 



1. INTRODUCTION 

The spectral energy distributions (SEDs) of BL Lac objects 
can be largely characterized by two broad peaks. The low 
and high energy peaks are commonly attributed to synchrotron 
emission and inverse Compton scattering, respectively. De- 
pending on the relative strengths of the two peaks, it is now 
customary to divide BL Lac objects into low frequency-peaked 
(LBL) and high frequency-peaked (HBL) sources (Giommi & 
Padovani 1994). Typically, HBLs are less variable, less polar- 
ized, and contain less dominant radio cores than LBLs (e.g., 
Laurent-Muehleisen et al. 1993; Perlman & Stocke 1993; Jan- 
nuzi, Smith & Elston 1994). However, the physical origin of the 
differences between the two classes of BL Lac objects remains 
a matter of debate. 

According to the unification scheme for BL Lac objects, 
LBLs and HBLs are observed at different viewing angles (see 
review by Urry & Padovani 1995). On the other hand, Sam- 
bruna, Maraschi, & Urry (1996) have shown that the typical 
multiwaveband SED of an HBL cannot be produced from the 
spectrum of an LBL simply by changing the viewing angle 
alone. They suggested that other physical parameters, such 
as the magnetic field, may be different in the two classes. 
New X-ray and radio surveys, such as those by Perlman et al. 

(1998) , Laurent-Muehleisen et al. (1998), and Caccianiga et al. 

(1999) , have shown that the distribution of peak frequencies in 
the SEDs of BL Lac objects is more uniform than previously 
thought. What is the main physical parameter (or set of param- 
eters) that control the SEDs of BL Lac objects? 

Recent observations of the host galaxies of BL Lac objects 
provide new clues to understanding the physical nature of these 



systems. Urry et al. (2000) and Scarpa et al. (2000) sys- 
tematically studied the morphologies of BL Lac objects using 
the Hubble Space Telescope (HST) and found that there is no 
significant difference between the host galaxies of HBLs and 
LBLs. Most of the hosts are normal, giant elliptical galaxies, 
with no obvious evidence for ongoing mergers or strong inter- 
actions with other galaxies. Urry et al. (2000) suggest that, if 
the mass of the central black hole (BH) in BL Lac objects cor- 
relates with the bulge luminosity of the host galaxy, as is the 
case in nearby, inactive galaxies (e.g., Magorrian et al. 1998; 
Kormendy & Gebhardt 2001), the Eddington ratios 5 in BL Lac 
objects must span a large range. A number of authors have sug- 
gested that the radio power in radio-loud quasars is controlled 
by the BH mass (e.g., Laor 2000; Lacy et al. 2001; McLure & 
Dunlop 2001; but see Ho 2002). If the host galaxies of BL Lac 
objects obey the BH mass/bulge luminosity relation, then the 
peak frequencies are evidently not controlled by the BH mass 
(Urry et al. 2000). This conclusion has been strengthened by 
the work of Falomo, Kotilainen & Treves (2002), Wu, Liu, & 
Zhang (2002), and Barth, Ho, & Sargent (2002a, b), who have 
obtained more robust BH mass estimates using the relation be- 
tween BH mass and the stellar velocity dispersion of the bulge 
(Gebhardt et al. 2000; Ferrarese & Merritt 2000; Tremaine et 
al. 2002). 

If BH mass is not the main parameter, what about the accre- 
tion rate? The importance of the role of accretion rates in BL 
Lac objects was originally recognized by Rees et al. (1982), 
who suggested that optically thin tori may power their central 
engines. A number of authors have recently revisited this is- 
sue in the context of BL Lac objects (Ghisellini & Celotti 2001; 
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The Eddington ratio refers to the ratio of the bolometric accretion luminosity to the Eddington luminosity, Ledd : 
mass in units of M@ . 
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FIG. 1 . — Distribution of peak frequencies for our sample of BL Lac objects. 
Note that there is no bimodality in the distribution. 

Maraschi 2001; Bottcher & Dermer 2002; Cavaliere & D'Elia 
2002; Cao 2002) and more generally radio-loud quasars (Lacy 
et al. 2001). Using a large sample of objects which covers a 
wide range of nuclear activity, from nearly inactive systems to 
classical Seyfert 1 nuclei and quasars, Ho (2002) shows that the 
conventional "radio-loudness" parameter is strongly inversely 
correlated with the mass accretion rate. 

This paper presents empirical evidence that the SEDs of BL 
Lac objects depend on the accretion rate. 

2. THE SAMPLE 

The sample of BL Lac objects is chosen according to the 
availability of information on (1) the luminosity of the host 
galaxy, (2) the luminosity of one or more emission lines, and 
(3) the fitted peak frequency and luminosity of the SED. The 
sample is limited mainly by the host galaxy observations. Rel- 
atively large ground-based optical imaging surveys of BL Lac 
objects have been performed by Abraham, McHardy, & Craw- 
ford (1991), Wurtz, Stocke, & Yee (1996), and Falomo & Koti- 
lainen (1999). More recently, an HST survey of 132 BL Lac 
objects was completed by Scarpa et al. (2000) and Urry et 
al. (2000). Whenever possible, we have given preference to 
the HST data because of the high resolution and the uniformity 
of that sample. Emission-line fluxes come primarily from the 
studies of Stickel, Fried, & Kiihr (1993), Rector et al. (2000), 
and Rector & Stocke (2001). Finally, we listed published val- 
ues of the peak frequencies and the corresponding peak lumi- 
nosities of the SEDs; these are obtained by fitting a logarithmic 
parabola to the multiwaveband, usually non-simultaneous, con- 
tinuum. 

The final sample, listed in Table 1, consists of 43 objects. 
For convenience, distance-dependent quantities assume Ho = 50 
km s" 1 Mpc" 1 and qo = 0. We note that the original classifica- 
tions by Urry et al. (2000), which used log^kev/^SGHz = -5.5 
as the criterion to separate the two types, contained some HBLs 
with low peak frequencies (i/ p k). Figure 1 shows the distribution 
of i/pk for our sample. There is no obvious bimodality in the dis- 
tribution; some objects have "intermediate" peak frequencies. 
For the subsequent analysis, we somewhat arbitrarily choose 
logi/pk = 15.1, roughly in the middle of the distribution, as the 
boundary between HBLs and LBLs. With this criterion, our 
sample has 17 HBLs and 26 LBLs. 
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NOTES: (a) Falomo & Kotilainen 1999; (b) Falomo & Ulrich 2000; (c) 
Pursimo et al. 2002; (d) Wolter et al. 1998; (e) Comastri, Molendi & Ghis- 
ellini 1995; (f) Landau et al. 1986. References for emission-line fluxes: (1) 
Rector et al. 2000; (2) Stickel, Fried, & Kiihr 1993. (3) Rector & Stocke 
2001; (4) Lawrence et al. 1996. The data for absolute magnitudes are taken 
from Urry et al. 2000, and the peak frequencies and luminosities are taken 
from Sambruna et al. 1996. The objects above and below the dividing are 
HBLs and LBLs, respectively. 

3. STATISTICAL PROPERTIES 

It is nontrivial to measure the thermal, accretion luminosity 
from BL Lac objects because the bulk of their emission comes 
from nonthermal processes associated with the jet. In this study, 
we adopt the strategy of using the observed luminosity of the 
emission lines to estimate the accretion luminosity of the disk. 
Emission lines themselves, however, are difficult to measure in 
BL Lac objects because they are both intrinsically weak and 
heavily diluted by the strong featureless continuum. Neverthe- 
less, they have been detected in a number of sources, and our 
sample was chosen with this application in mind. 

In our sample, the most commonly detected emission lines 
are Mgn A2800, [On] A3727, H/3, [Our] A5007, and Ha. To 
estimate the total luminosity due to emission lines, we follow 
the method of Celotti, Padovani & Ghisellini (1997), who com- 
pare the luminosity of the observed lines with the fractional 
contribution these lines make to the total line luminosity, as 
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FIG. 2. — The distribution of accretion rates in BL Lac objects. The dotted 
histogram represent the HBLs and solid histograms show the LBLs. 

determined from the relative line ratios computed by Fran- 
cis et al. (1991) for a composite quasar spectrum. As in 
Celotti et al. (1997), we assume that the total line luminosity 
(^lines) = 556(LLya)- The line luminosities are given in Table 1. 

With the assumption that the bulk of the line luminosity is 
photoionized by the accretion disk (Netzer 1990), Lii nes should 
be proportional to the total disk luminosity. We define a "line 
accretion rate" and its dimensionless form as follows: 



^ _ ^lines t 



A = 



^Edd 



(1) 



The relation between A and the dimensionless accretion rate m 
is given by equation (2). We determine BH masses from the 
absolute /?-band magnitudes of the host galaxies, using the re- 
lation determined by McLure & Dunlop (2001) from a study 
of active galaxies: logm BH = -0.5M R -2.96. This relation is 
quite similar, but not identical, to that given by Kormendy & 
Gebhardt (2001) for inactive galaxies. 

The above estimates of the accretion rates may carry signif- 
icant uncertainties that are difficult to quantify. While the for- 
mal errors on the luminosities of individual emission lines are 
usually quite modest, what is less certain is the conversion to 
total line luminosity, as this step implicitly assumes that the 
relative line strengths in BL Lac objects are the same as in typ- 
ical quasars. Thus, our values of Lii nes , and hence A and A, are 
potentially affected by systematic errors. The values of A are 
additionally affected by the relatively large scatter in the BH 
mass versus bulge luminosity relation (^0.59 dex; McLure & 
Dunlop 2001). Despite the potentially large uncertainties on 
any individual measurement of A or A, and their inferred accre- 
tion rates, we would like to emphasize that the relative values 
of these quantities, such as those deduced for HBLs and LBLs 
discussed below, should be more robust. 

3.1. Accretion Rates 

Figure 2 shows that the accretion rates vary over a rather 
wide range in LBLs, whereas they are more restricted in HBLs. 
The mean values of the accretion rates are (log A) = 22.70 ± 
0.91(rms dispersion) and (log A) = -3.82 ±0.68 in LBLs and 
(log A) = 21.69 ±0.56 and (log A) = -4.48 ±0.62 in HBLs. 
These results directly show that the accretion rate may control 
the SEDs of BL Lac objects. It is important to note that the dis- 
tributions of A and A significantly overlap between HBLs and 
LBLs. This clearly reflects the continuous distribution of peak 
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FIG. 3. — The correlations among the parameters u p ^, Lp^, and A. The solid 
and open symbols denote LBLs and HBLs, respectively. 

frequencies in Figure 1 . 

With the assumption that the line luminosity traces the disk 
luminosity, the distribution of A implies that the vast majority 
of BL Lac objects are highly sub-Eddington systems. In such 
a regime, accretion is thought to proceed via an optically thin 
"advection-dominated" accretion flow (ADAF; see Narayan, 
Mahadevan, & Quataert 1998). From Mahadevan (1997), the 
total luminosity from an ADAF is oc a~ 2 m BH m 2 , where 
m = M/M Edd , M Edd = 1.39 x I0 is m m /'q-i (g s _1 ), 77 = O.lr/-! is 
the accretion efficiency, and a is the canonical viscosity param- 
eter. For a reprocessing efficiency £, Lii nes = £Ldisk and 



m = 2.17 x 10" z ao.3?_i A 



1/2,1/2 



-4 > 



(2) 



where cto.3 = a/0.3, £-1 = £/0.1 and A-4 = A/10" 4 . We assume 
that £ is of the same order as the covering factor of the line- 
emitting clouds (^0.1) since the total energy emitted by the 
clouds is simply the energy absorbed by them (Netzer 1990). 

An optically thin ADAF can only exist below a critical value 
of the accretion rate m c « a 2 (Narayan et al. 1998). For a plau- 
sible value of a = 0.3, m c w 0.1. We find (logm) =-1.5 ±0.34 
for LBLs and (logm) = -2.0 ± 0.31 for HBLs. This indicates 
that BL Lac objects are in the ADAF regime, and that LBLs 
typically have significantly higher accretion rates than HBLs. 
These results are in qualitative agreement with the findings of 
Ghisellini & Celotti (2001) and Cavaliere & D'Elia (2002). We 
note that equation (2) may not be applicable to some of the 
LBLs with the highest values of A because they have m ~ m c ; 
these objects straddle the (uncertain) transition from a standard 
disk to an ADAF (cf. Cao 2002). 
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A for HBLs and LBLs can be understood in terms of different 
radiative efficiencies in the jet. 

Let us assume that the energy density of the random com- 
ponent of the magnetic field is proportional to the energy den- 
sity of the relativistic electrons in the jet. Then, B 2 oc (E e )n e , 
where (E e ) is the mean energy and n e is the number density 
of the relativistic electrons. With the help of Mj et ~ R\n e m p c, 
where Rq is the radius of the jet and m p is the proton mass, 

we have B oc M-^R^ 1 and n e oc M^R^. If the relativistic elec- 
trons have a power-law distribution of energies N = NoE' 1 , the 
emissivity at the peak frequency i/ k in the jet comoving frame 

is given by e pk oc NqB^^u'^' 2 oc n e B<T +1 >/ 2 z^ 7>/2 (Pachol- 

czyk 1970), and it follows that e pk oc Rf^ 2 ^ 5 ^ 4 ^' 1 . 
The radiative luminosity in the jet comoving frame will then be 
given by L' ai oc R\ (, pk£pk ) cx R^ 2 ^ 4 ^- The ra- 
tio of the radiative to the kinetic luminosity will be proportional 
to r -1 ^^ 7 ^ 4 ^ ~'" 2 , which suggests that the higher peak fre- 
quencies lead to the higher radiative efficiencies. 

For HBLs, we consider the limiting case that most of the 
energy is radiated. By setting L' rad ~ L B z 
M~(7+i)/ 2 (3-7) ) an( j the observed peak frequency is 



we have i/' pk oc 



= Vv' pk oc M-(7+D/2(3-7) oc A-^ +1 >/ 4 ^. 



(7) 



3.2. Accretion Rates and SEDs 

Figure 3 shows the dependence of the peak frequency O p k) 
and peak luminosity (L pk ) on A. We find that L ?k correlates 
strongly with A; the Pearson's correlation coefficient is r = 0.76, 
significant at a level of >99.999%. Treating A as the indepen- 
dent variable, a linear regression fit gives 

logL pk (erg s" 1 ) = (30.88 ± 1.99) + (0.64 ±0.09) log A. (3) 

An obvious interpretation of the L p k- A correlation is that the jet 
power is directly coupled to the mass accretion rate of the disk. 
The disk-jet "connection" in active nuclei has been discussed in 
a number of contexts (e.g., Rawlings & Saunders 1991; Falcke 
& Biermann 1995; Xu, Livio & Baum 1999; Ho & Peng 2001), 
and the present correlation provides additional support for this 
picture because the peak luminosity of the SED of a BL Lac 
object represents most of the radiated energy from the jet. 

The relation between v vk and A is more complex. The 
HBL group and the LBL group each seems to define its 
own sequence. After excluding three strongly deviant points 
(0922+749, 1229+643 and 1458+224) 6 , HBLs appear to show 
a moderately strong correlation with r = -0.6 at a significance 
of 98%. The best-fit linear regression gives 

log j/p k (Hz) = (81 .65 ± 25 .61) - (3.02 ± 1 . 19) log A. (4) 

The peak frequency formally scales with the accretion rate as 
f p k oc A" 3 oc M~ 6 , but we note that this is somewhat uncertain 
because of the limited dynamic range in A. For LBLs, we have 

log fp k (Hz) = (19.66 ± 2.55) - (0.25 ± 0. 1 l)log A, (5) 
with r = -0.42 and a significance of 97%. This implies z/ p t oc 
A -0 25 oc M"° 5 , much shallower than in HBLs. Bearing in mind 
the above caveat on the HBLs, the different slopes of the v v \- A 
relation may indicate differences in the physical processes as- 
sociated with the jet, as discussed in the next section. It is in- 
teresting to note that HBLs and LBLs also delineate different 
loci, separated by a sharp break, when is plotted against the 
radio luminosity (Costamante et al. 2001). 

4. DISCUSSION 

According to Xu et al. (1999), the accretion rate through the 
disk (M) is related to the mass flux into the jet (Mj et ) by 

M iet /M~(H/R) 2S -\ (6) 

where H is the height of the disk at radius R and S is a con- 
stant from 1.7-3.4. From the distribution of A, we inferred 
that both HBLs and LBLs contain ADAFs. Since the struc- 
ture of ADAFs is quasi-spherical, H/R ~ 1 or Mj et ~ M. As 
Blandford & Begelman (1999) argued, ADAFs inevitably lead 
to outflows. If energy is extracted from the spin of the BH via 
the Blandford-Znajek (1977) mechanism, Lbz = £bzMc 2 , where 
€bz = 1 .8 x 10~ 2 and the BH spin is assumed to be maximal (Ar- 
mitage & Natarajan 1999). The BZ process has been strength- 
ened by a recent observation of Fe Ka line profile in MCG 
6-30-15 (Wilms et al. 2001). This energy will be channeled 
into two parts: the kinetic luminosity Lyn = Mj et Tc 2 /2 (T is the 
Lorentz factor of the jet) and the radiative luminosity L' md (seen 
in the comoving frame of the jet; e.g., Sikora et al. 1997). Here, 
we explore the possibility that the relative partition of the spin 
energy determines the shape of the SED. In particular, we in- 
vestigate whether the different functional dependence of v v ^ on 

6 The emission-line equivalent widths of 0922+749 and 1229+643 formally exceed the (admittedly arbitrary) limit of 5 A. The value of u pk for 0922+7 '49 was obtained 
from fitting only three points (Sambrana et al. 1996) and is therefore quite uncertain. There is no obvious explanation for 1458+224. 



Here, we have assumed that the Doppler factor V does not vary 
strongly with M. We see that v p k depends sensitively on A for 
a given 7. For a typical value of 7 « 2.6-2.7 (Kirk, Melrose 
& Priest 1994), we get z/ pk oc A" 2 3 to A" 31 . This is consistent 
with the empirical results for HBLs (eq. 4). Despite the appar- 
ent agreement, however, we must regard this result with some 
caution because of the extreme requirement that the bulk of the 
jet power emerge as radiative energy. Celotti & Fabian (1993) 
find that most jet sources do not radiate efficiently, although 
there are some objects with high radiative efficiencies (see their 
Fig. 2). 

Since LBLs have lower peak frequencies, their radiative ef- 
ficiencies are lower, and we postulate that in this case the ki- 
netic luminosity dominates over the radiative luminosity. The 
maximum energy of the relativistic electrons will be given by 
the balance between energy loss and gain. The acceleration 
timescale is r acc ~ l{E e )/c oc B~ 1 E", where l(E e ) = IqE" is the 
mean free path and n = 1 in the Bohm limit (Inoue & Takahara 
1996). The timescale for energy loss due to inverse Compton 
scattering of external photons is ti oss oc (M e xt^e) _1 > where the en- 
ergy density of external photons u ext oc Ldi.sk (Sikora, Begelman 

& Rees 1994). We thus have E e ^ 
peak frequency is 



B/Ldisk) 1//( " +1) , and the 



u pk oc VBE 2 e umx ocM ( "" 5) / 2( " +1) oc A ( "- 5) / 4( " +1 \ 



(8) 



where we used Ldisk oc M 2 . For the Bohm limit in = 1) or 
not very efficient acceleration (e.g., n = 2), the relation will 
be z/pk oc A" 0,25 to A" 0,50 . This, too, is roughly consistent with 
the empirical correlation for LBLs (eq. 5). In this scenario, 
which is similar to that discussed by Ghisellini et al. (1998) 
and Bottcher & Dermer (2002), the peak frequencies in LBLs 
are determined by cooling of relativistic electrons due to inverse 
Compton scattering of photons from the disk or reflected off of 
the broad-line emitting clouds. 
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Finally, we would like to note that our analysis hinges on 
the assumption that the emission lines in BL Lac objects are 
photoionized by the disk, such that line luminosity reliably 
traces the accretion luminosity. Based on spectrophotometric 
monitoring of a sample of high-redshift, high-luminosity active 
galactic nuclei, Perez, Penston & Moles (1989) find that the 
emission lines vary on timescales shorter than expected from 
photoionization models. They suggest that the broad-line re- 
gion in these objects may be anisotropic. However, the vari- 
ability behavior of the broad optical line Ha in BL Lac itself 
appears to support that "Ha emission could be powered by 
thermal radiation from an accretion disk without significantly 
affecting the shape ot polarization of the optical continuum" 
(Corbett et al. 1996, 2000). If it turns out that a significant 
fraction of the line luminosity is powered by the jet instead of 
the disk, then clearly the interpretation of Figure 3 needs to be 
reevaluated. However, the conclusion that BL Lac objects are 
highly sub-Eddington, advection-dominated systems would not 
change; in fact, it would only be strengthened because the in- 
ferred accretion luminosity would be even lower than we have 
assumed. 

5. CONCLUSIONS 

Using a sample of BL Lac objects with information on their 
host galaxies and emission-line luminosities, we show that 
LBLs have significantly higher accretion rates and Eddington 
ratios than HBLs. Both classes are highly sub-Eddington sys- 



tems that may be accreting via an ADAF. We find that the peak 
luminosity of the SED correlates significantly with the accre- 
tion rate, lending strong support to the idea that the jet and ac- 
cretion disk are closely coupled. 

Furthermore, we present evidence that the peak frequency of 
the SED also correlates with the accretion rate, although the 
functional dependence between the two parameters appears to 
be quite different between HBLs and LBLs. We argue that these 
empirical trends can be qualitatively explained by invoking dif- 
ferent radiative efficiencies in the jets of the two classes of ob- 
jects. HBLs, which have lower accretion rates, evidently man- 
age to convert more of their jet kinetic power into radiation. By 
contrast, in LBLs, which have higher accretion rates, a greater 
fraction of the jet power remains in kinetic form. However, the 
underlying causal connection between the radiative efficiency 
of the jet and the accretion rate, and the manner in which en- 
ergy is channeled into relativistic electrons (Ghisellini, Celotti 
& Costamante 2002), remain unclear. 

It would be of interest to perform spectrophotometric moni- 
toring of nearby BL Lac objects in order to establish more con- 
clusively the origin of their broad emission lines. 
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